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ABSTRACT 

Measurements of ultraviolet line fluxes from Space Telescope Imaging Spectrograph and Far- 
Ultraviolet Spectroscopic Explorer spectra of the K2-dwarf e Eri are reported. These are used 
to develop new emission measure distributions and semi-empirical atmospheric models for 
the chromosphere and lower transition region of the star. These models are the most detailed 
constructed to date for a main-sequence star other than the Sun. New ionisation balance cal- 
culations, which account for the effect of finite density on dielectronic recombination rates, 
are presented for carbon, nitrogen, oxygen and silicon. The results of these calculations are 
significantly different from the standard Arnaud & Rothenflug ion balance, particularly for 
alkali-like ions. The new atmospheric models are used to place constraints on possible First 
Ionisation Potential (FIP) related abundance variations in the lower atmosphere and to discuss 
limitations of single-component models for the interpretation of certain optically thick line 
fluxes. 
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1 INTRODUCTION 

Although it has long been known that cool main-sequence stars 
have hot outer atmospheres the question of how these regions 
are heated remains unanswered. This is, in part, due to the com- 
plexity of stellar atmospheres arising from the interplay of sub- 
photospheric convection, rotation and magnetic fields - in order to 
distinguish between different physical processes quantitative anal- 
yses and models are required. 

Atmospheric models fall into two broad categories, those 
which are constructed from theory (e.g. those developed by 
Fontenla, Avrett & Loeser 1993, 2002 or Cuntz et al. 1999) and 
those which are constructed semi-empirically, based on a combi- 
nation of observations and a few elementary physical assumptions 
(e.g. the Vernazza, Avrett & Loeser 1973, 1976, 1981 solar mod- 
els). Semi-empirical models are vital since it is only by comparison 
with properties deduced from observations that theoretical models 
can be tested. 

In this paper new semi-empirical models are developed for the 
chromosphere and lower transition region of the K2-dwarf e Eri. 
This complements our recent studies of the outer atmosphere in 
which the mean electron pressure was measured (Jordan et al. 
2001a); line widths were measured and used to investigate the non- 
thermal velocity fields within the atmosphere and place constraints 
upon the possibility of coronal heating by wave dissipation (Sim & 
Jordan 2003a); and emission measure distributions were used to de- 
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rive new information about the area occupied by emitting material 
in the upper transition region (Sim & Jordan 2003b). 

Previously, semi-empirical models have been constructed for 
e Eri by Kelch (1978), Simon, Kelch & Linsky (1980) and Thatcher, 
Robinson & Rees (1991). The Kelch (1978) model was based on 
observations of the Ca II and Mg II resonance lines and described 
the upper photosphere and chromosphere (extending to an electron 
temperature logT e [K] = 4.3). Simon et al. (1980) used ultraviolet 
observations of Si II, Si III and C II recorded with the International 
Ultraviolet Explorer (IUE) to develop the upper parts of the Kelch 
(1978) model. Their model extends to logT e [K] = 4.7. Thatcher 
et al. (1991) used infrared lines of Ca II, the Na D doublet and lines 
from the H B aimer series to further develop the chromospheric part 
of the Kelch (1978) model and extended the model to log T e [K] = 
5.6. 

The new models presented here are based on ultraviolet ob- 
servations made with the Space Telescope Imaging Spectrograph 
(STIS) on board the Hubble Space Telescope and with the Far- 
Ultraviolet Spectroscopic Explorer (FUSE) satellite (see Section 
2). The high signal-to-noise ratio, high spectral resolution and large 
combined spectral coverage of these data are ideal for quantitative 
modelling; they provide much tighter constraints on the model than 
did the IUE data and allow much more robust conclusions to be 
drawn. In particular, the improved measurements of the transition 
region pressure, together with radiative transfer calculations, lead 
to more reliable models from the photosphere to the transition re- 
gion. 

An important limitation of the available spectroscopic data is 
that, in common with all such observations of main-sequence stars 
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other than the Sun, they contain no direct information on the spa- 
tial inhomogeneities (because the source is unresolved). Although 
we have been able to deduce information about the inhomogene- 
ity in the upper transition region (Sim & Jordan 2003b), the same 
method cannot yet be applied to the deeper parts of the atmosphere 
since the local means of energy deposition is not understood. Here, 
for the most part, the simplest assumption, that of homogeneity, is 
made - thus the models are representative of the "average" atmo- 
sphere. Although the interpretation of "average" models requires 
considerable care, they are extremely important since they con- 
tain the most information that can be extracted directly from the 
data alone and theoretical models must be consistent with the mean 
properties. Two-component theoretical models of a K2 V star chro- 
mosphere have been made by Cuntz et al. (1999), using empirical 
relations between the area covered by magnetic fields and stellar 
rotation periods. We also derive a two-component by applying the 
emitting area that we found at 10 5 K (Sim & Jordan 2003b) as an 
upper limit to regions at lower temperatures - solar observations 
show that the area occupied by supergranulation cells remains con- 
stant between ~ 2 x 10 4 K and 10 5 K, before decreasing in the 
chromosphere. 

Another limitation to the modelling is the accuracy of the 
atomic data used - for many of the measured line fluxes, the ob- 
servational uncertainty is small compared to possible errors in the 
atomic rates. In particular, the standard sets of ionisation balance 
calculations used for transition region lines are inadequate for mod- 
elling to the precision of the STIS data. In this paper a new ionisa- 
tion balance calculation is presented in which rates that account for 
variations in the dielectronic recombination rate with density are 
used; although not definitive, this calculation highlights the need 
for further work on recombination rates at non-zero density. 

The ultraviolet observations and measurements of the emis- 
sion line fluxes are discussed in Section 2. In Section 3 the new 
ionisation balance calculations are discussed. They are used in the 
construction of emission measure distributions (EMDs) in Section 
4. A preliminary single-component atmospheric model is devel- 
oped in Section 5 and improvements to this model are discussed 
in Section 6. In Section 7, the simple two-component model is pre- 
sented and conclusions are drawn in Section 8. 



2 OBSERVATIONS 
2.1 Flux measurements 

The models developed are based on measurements of fluxes for a 
variety of emission lines in STIS and FUSE spectra of e Eri. Ta- 
ble 1 gives the measured fluxes for lines in the STIS spectra (see 
Jordan et al. 2001b for details of our STIS spectra). These were 
obtained by direct integration across line profiles and, where possi- 
ble, by performing single-Gaussian fits. (The line widths and shifts 
deduced from Gaussian fits have been published elsewhere: Sim & 
Jordan 2003a.) The errors given are tolerances based on the signal- 
to-noise, the uncertainty in the background level and, where appli- 
cable, the goodness of the Gaussian fits. Four of the lines listed in 
Table 1 are affected by blending. Three of these (O IV 1404.8 A, 
Si II 2335.3 A and S IV 1404.8 A) have been discussed by Jor- 
dan et al. (2001a). The other, Si II 1260.4 A, is suspected to be a 
close blend with C I 1260.6 A. The Al III 2 Si /2 - 2 Pi/ 2 line at 
1862.8 A it too weak to measure. This suggests that the Al III 
1854.7-A line is also blended, probably with an iron line. (In the 
optically thin limit the Al III 1854.7-A line should have double the 



flux of the 1862.8-A line, but the observed 1854.7-A feature is more 
than twice as strong as that at 1862.8 A.) 

It is noteworthy that when the IUE fluxes reported for e Eri by 
Jordan et al. (1987) are compared with the new fluxes from STIS 
there is no evidence of significant differences. Since the IUE and 
STIS data were recorded seventeen years apart this argues against 
significant long term variability in the atmosphere. 

The STIS fluxes are supplemented by measurements from 
FUSE spectra of e Eri. The FUSE observations are discussed by 
Redfield et al. (2002). The lines in the FUSE spectrum that are used 
here, and our measurements of their fluxes, are listed in Table 2. As 
for the STIS measurements, the fluxes were obtained by a com- 
bination of direct integration and single-Gaussian fits. The FUSE 
instrument has four channels which record spectra, 1-LiF, 2-LiF, 1- 
SiC and 2-Sic. Most of the FUSE line fluxes have been measured 
from spectra recorded with the 1-LiF channel (this channel is used 
to guide the FUSE line-of-sight). Unfortunately, the S VI 933.4- A, 
C III 977.0-A and H I 972.5-A lines do not fall in the 1-LiF spec- 
tra. Therefore, the fluxes for these lines have been measured from 
1-SiC spectra which are susceptible to losing photons due to chan- 
nel drift (Sankrit et al. 2001). Measurements of the strong O VI line 
fluxes (which occur in both the 1-LiF and 1-SiC spectra) were used 
to correct the 1-SiC line fluxes for channel drift but the measure- 
ment uncertainties are still largest for the 1-SiC measurements. 

Redfield et al. (2002) have published independent measure- 
ments of these line fluxes in the e Eri spectrum. In almost every 
case, the fluxes here agree with those given by Redfield et al. (2002) 
to within the quoted errors. The only exception is the C III 1175- 
A multiplet; in this case, Redfield et al. (2002) report a flux which 
is smaller by ~ 20 per cent than that found here, while the errors 
are only ~ 10 per cent. It is noted that their fit (their fig. 10) does 
appear to underestimate the peak fluxes in this multiplet. However, 
this difference is not significant given the other sources of uncer- 
tainty in the modelling. The total C III 1 175-A flux measured from 
STIS (Table 1) and that derived from FUSE (Table 2) are very sim- 
ilar (they agree to within the measurement uncertainties). This ar- 
gues against short term variations in the emission line fluxes since 
the STIS and FUSE data were recorded nine months apart. 



2.2 Interstellar absorption 

Interstellar absorption can be clearly identified in the line profiles 
of several of the resonance lines. This is important since it allows 
information to be derived about the local interstellar medium (see, 
for example, Wood et al. 2000 or Redfield & Linsky 2002 for such 
studies involving spectra of e Eri) but also because it must be taken 
into consideration when analysing the spectrum of the star. The 
largest effect is in lines of the H I Lyman series (this H I absorption 
has been studied by Dring et al. (1997) who used it to determine 
the interstellar H I and D I column densities). However, there is 
also significant absorption in the Mg II h and k lines and several 
transitions involving the ground states of C II and Si II. The C II 
and Si II lines affected are listed in Table 3 which also gives esti- 
mated line centre opacities (to) for interstellar absorption in each 
line. These have been computed with the C/H abundance for the 
interstellar medium given by Cardelli et al. (1996) (gas phase C/H 
ratio of 1.4 x 10~ 4 ) and assuming a solar photospheric value for the 
relative Si/C abundance (Grevesse & Sauval 1998). The interstel- 
lar H I column density was taken from Dring et al. (1997) and the 
interstellar line width parameter was estimated from the observed 
interstellar absorption in the Mg II lines. For the calculations of 
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Table 1. Emission line fluxes measured from the e Eri STIS spectra. Ao is the rest wavelength, "bl." indicates that the line is part of a blend: in such cases the 
total measured flux in the blend is given as an upper limit to the line flux. The rest wavelengths are taken from Griesmann & Kling (2000) [C IV and Si IV], 
Kaufman & Martin (1993) [S IV] and Kurucz & Bell (1995) [all other lines]. 



Ion 


Transition 


Ao 


Flux at Earth 








(A) 


/-in— IS — 2 — 1 \ 

(10 J ergs cm z s 1 ) 


Hi 


ls 2 S 


-2p 2 P 


1215.671 


160 ± 1 


CII 


2s 2 2p 2 P 1/2 - 


-2s2p 2 2 D 3/2 


1334.532 


1.88 ± 0.06 


CII 


2s 2 2p 2 P 3/2 -2s2p 2 2 D 5/2j3/2 


1335.68 bl. 


4.84 ±0.10 


C III 


2s2p 3 Pi 


- 2p 2 3 P 2 


1174.933 


0.41 ± 0.06 


C III 


2s2p 3 P 


-2p 23 P! 


1175.263 


0.41 ± 0.06 


C III 


2s2p 3 P 2 ,i 


-2p 2 3 P 2 .! 


1 175.69 bl. 


1.49 ± 0.06 


C III 


2s2p 3 Pi 


-2p 2 3 P 


1175.987 


0.41 ± 0.06 


C III 


2s2p 3 P 2 


-2p 2 3 Pi 


1176.370 


0.46 ± 0.06 


C IV 


2s 2 S 1/2 - 


-2p 2 P 3/2 


1548.204 


5.59 ± 0.11 


C IV 


2s 2 Si /2 - 


-2p 2 Pi /2 


1550.781 


2.82 ± 0.07 


NIV 


2s 2 !S - 


•2s2p 3 Pi 


1486.496 


0.045 ± 0.015 


N V 


2s 2 S 1/2 - 


-2p 2 P 3/2 


1238.821 


0.94 ±0.03 


N V 


2s 2 S 1/2 - 


-2p 2 Pi/ 2 


1242.804 


0.46 ± 0.02 


O III 


2s 2 2p 2 3 Pi 
2s 2 2p 2 3 P 2 


-2s2p 3 5 S 2 


1660.809 


0.058 ± 0.015 


O III 


-2s2p 3 5 S 2 


1666.150 


0.121 ± 0.015 


OIV 


2s 2 2p 2 P 1/2 


-2s2p 2 4 P 1/2 


1399.780 


0.054 ± 0.005 


OIV 


2s 2 2p 2 P 3/2 


-2s2p 2 4 P 5/2 


1401.157 


0.181 ± 0.005 


OIV 


2s 2 2p 2 P 3/2 


-2s2p 2 4 P 3/2 


1404.806 


0.069 ± 0.005 bl. 


OIV 


2s 2 2p 2 P 3/2 


-2s2p 2 4 P 1/2 


1407.382 


0.052 ± 0.005 


O V 


2s 2 1 S - 


2s2p 3 Pj 


1218.344 


0.68 ± 0.04 


O V 


2s2p x Pi 


-2p 21 D 2 


1371.296 


0.105 ± 0.005 


Mg II 


3s 2 S 1/2 - 


-3p 2 P 3/2 


2796.352 


305 ±6 


Mg II 


3s 2 S 1/2 - 


-3p 2 Pi/ 2 


2803.531 


232 ±5 


Aim 


3s 2 S 1/2 - 


-3p 2 P 3 /2 


1854.716 


0.51 ±0.06 bl? 


Sill 


3s 2 3p 2 P 1/2 


-3s3p 2 2 P 3/2 


1190.416 


0.10 ±0.03 


Sill 


3s 2 3p 2 P 1/2 


-3s3p 2 2 P 1/2 


1193.290 


0.09 ± 0.04 


Sill 


3s 2 3p 2 P 3/2 


-3s3p 2 2 P 3/2 


1194.500 


0.18 ± 0.03 


Sill 


3s 2 3p 2 P 3/2 


-3s3p 2 2 P 1/2 


1197.394 


0.10 ±0.04 


Sill 


3p 2 Pi/ 2 - 


-3d 2 D 3/2 


1260.422 


0.191 ±0.03 bl. 


Sill 


3p 2 P 3/2 ■ 


-3d 2 D 5/2 


1264.738 


0.505 ± 0.05 


Sill 


3p 2 P3/ 2 - 


-3d 2 D 3/2 


1265.002 


0.172 ± 0.03 


Sill 


3s 2 3p 2 P 1/2 


-3s3p 2 2 S 1/2 


1304.370 


0.133 ± 0.008 


Sill 


3s 2 3p 2 P 3/2 


-3s3p 2 2 S 1/2 


1309.276 


0.222 ± 0.010 


Sill 


3p 2 Pi/ 2 


-4s 2 S 1/2 


1526.707 


0.508 ± 0.005 


Sill 


3p 2 P 3 /2 


-4s 2 Si /2 


1533.431 


0.607 ±0.005 


Sill 


3s 2 3p 2 P 1/2 - 


-3s3p 2 2 D 3/2 


1808.013 


6.39 ± 0.06 


Si II 


3s 2 3p 2 P 3/2 - 


-3s3p 2 2 D 5/2 


1816.928 


10.42 ± 0.07 


Si II 


3s 2 3p 2 P 3/2 - 


-3s3p 2 2 D 3/2 


1817.451 


3.20 ±0.06 


Sill 


3s 2 3p 2 P 1/2 


-3s3p 2 4 P 1/2 


2335.123 


3.31 ± 0.07 


Sill 


3s 2 3 P 2 P 3/2 


-3s3p 2 4 P 5/2 


2335.321 


2.23 ± 0.07 bl. 


Sill 


3s 2 3p 2 P 3/2 


-3s3p 2 4 P 1/2 


2350.891 


3.04 ±0.07 


Sim 


3s 2 1 S - 


■3s3p x Pi 


1206.500 


3.60 ±0.10 


Sim 


3s3p 3 Pi 


-3 P 2 3 P 2 


1294.545 


0.052 ± 0.005 


Sim 


3s3p 3 P 


- 3p 2 3 Pi 


1296.726 


0.041 ± 0.005 


Sim 


3s3p 3 P 2 ,i 


-3p 2 3 P 2 ,i 


1298.93 bl. 


0.167 ±0.005 


Sim 


3s3p 3 Pi 


-3 P 2 3 P 


1301.149 


0.028 ± 0.005 


Sim 


3s3p 3 P 2 


~3p 2 3 Pi 


1303.323 


0.062 ± 0.005 


Sim 


3s 2 i So - 


-3s3p 3 Pi 


1892.030 


1.25 ± 0.05 


Si IV 


3s 2 S 1/2 - 


-3p 2 P 3/2 


1393.760 


2.30 ±0.06 


Si IV 


3s 2 S 1/2 - 


-3p 2 Pi/ 2 


1402.773 


1.21 ± 0.05 


SII 


3s 2 3p 34 S 3/2 


-3s3p 4 4 P 1/2 


1250.584 


0.029 ± 0.010 


Sn 


3s 2 3p 3 4 S 3/2 


-3s3p 4 4 P 3/2 


1253.811 


0.066 ± 0.005 


SII 


3s 2 3p 34 S 3/2 


-3s3p 4 4 P 5/2 


1259.519 


0.099 ± 0.005 


S IV 


3s 2 3p 2 P 1/2 


-3s3p 2 4 P 1/2 


1404.808 


0.069 ± 0.005 bl. 


S IV 


3s 2 3p 2 P 3/2 


-3s3p 2 4 P 5/2 


1406.016 


0.032 ± 0.007 


SIV 


3s 2 3p 2 P 3/2 


-3s3p 2 4 P 3/2 


1416.887 


0.042 ± 0.007 


CI I 


3p 5 2 Pi /2 -3p 4 ( 3 P)4s 2 P 1/2 


1351.656 


0.198 ± 0.005 
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Table 2. Fluxes for emission lines in the e Eri FUSE spectra. Unless otherwise indicated, all lines were measured in the 1-LiF spectra. Ao is the rest wavelength. 
The rest wavelengths are taken from Kurucz & Bell (1995). 



Ion 


Transition 




Flux at Earth 






(A) 


/-in— 1 ^ —9 —Is 

(10 1J ergs cm A s 1 ) 


Hi" 


ls 2 S-4p 2 P 


972.537 


0.63 ±0.08 


Hi 


ls 2 S-3p 2 P 


1025.722 


3.66 ±0.12 


Cn 


2s2p 2 4 P 1/2 -2p 3 4 S 3/2 


1009.858 


0.012 ±0.007 


C II 


2s2p 2 4 Po /2 -2p 3 4 So /2 


1010.083 


0.027 ± 0.007 


C II 


2s2p 2 4 P=,, -2p 3 4 So /9 

r O / Z C O j Z 


1010.371 


0.046 ± 0.007 


C II 


2s 2 2d 2 Pi to - 2s2d 2 2 Si /o 


1036.337 


0.120 ± 0.014 


C II 


2s 2 2p 2 P 3 / 9 -2s2p 2 2 S, io 

r o j Z r 1 / Z 


1037.018 


0.291 ± 0.008 


C III" 


2s 2 1 S -2s2p x Pi 


977.020 


4.67 ± 0.10 


C III 


2s2p 3 Pi - 2p 2 3 P 2 


1174.933 


0.435 ±0.015 


C III 


2s2p 3 P - 2p 2 3 Pi 


1175.263 


0.425 ±0.015 


C III 


2s2p 3 P 2 , 1 -2p 2 3 P 2>1 


1175.63 bl. 


1.519 ±0.015 


C III 


2s2p 3 Pi - 2p 2 3 P 


1175.987 


0.405 ±0.015 


C III 


2s2p 3 P 2 - 2p 2 3 Pi 


1176.370 


0.467 ±0.015 


Nm 


2s 2 2p 2 P 3/2 - 2s2p 2 2 D 3/2j5/2 


991.54 bl. 


0.33 ±0.08 


O VI 


2s 2 S 1/2 -2p 2 P 3/2 


1031.912 


4.59 ±0.16 


O VI 


2s 2 S 1/2 -2p 2 P 1/2 


1037.613 


2.26 ±0.08 


Sim 


3s3p 3 P - 3s3d 3 Di 


1108.358 


0.020 ±0.020 


Sim 


3s3p 3 Pj -3s3d 3 D lj2 


1 109.96 bl. 


0.065 ± 0.020 


Sim 


3s3p 3 P 2 -3s3d 3 Di, 2i3 


1113.20 bl. 


0.121 ±0.020 


S IV 


3s 2 3p 2 Pi /2 -3s3p 2 2 D 3/2 


1062.678 


0.038 ±0.020 


S IV 


3s 2 3p 2 P 3/2 -3s3p 2 2 D 5/2 


1072.996 


0.056 ±0.020 


S VI a 


3s 2 S 1/2 -3p 2 P 3/2 


933.378 


0.24 ±0.08 



a Measured from 1-SiC spectra. 



Table 3. Lines of C II and Si II affected by interstellar absorption and their 
calculated interstellar line centre opacities (to) The absolute values of to 
are uncertain but the relative values are reliable (see text). 



Ion 


Ao (A) 


TO 


Cn 


1334.5 


5.6 


CII 


1036.3 


4.4 


Si II 


1 190.4 


1.7 


Si II 


1193.3 


3.5 


Si II 


1260.4 


7.5 


Si II 


1304.4 


0.2 


Sill 


1526.7 


1.0 


Sill 


1808.0 


2 x 10~ 2 



to, it has been assumed that all the interstellar C and Si is singly 
ionised - thus the listed values are upper limits. 

The relative values of To suggest that the 1260.4- A line should 
show the greatest interstellar absorption of the Si II lines and that 
there should be very little effect on the 1526.7, 1304.4 or 1808.0- 
A lines. This is consistent with the observed line profiles - of the 
Si II lines, only the 1260.4- A transition shows evidence of self- 
reversal (in addition to line blending). The C II resonance lines 
at 1334.5 A (observed with STIS) and 1036.3 A (observed with 
FUSIi) are expected to have similar optical depths due to interstel- 
lar absorption and both show evidence of absorption. 



3 IONISATION BALANCE CALCULATIONS 

Ionization equilibrium populations are needed for quantitative anal- 
yses of stellar atmospheres. Above the chromosphere, for most ions 
the dominant ionisation processes are collisional; either direct col- 
lisional ionisation or collisional excitation followed by autoionisa- 



tion (Goldberg, Dupree & Allen 1965). Photoionisation is usually 
not important. Both radiative recombination and dielectronic re- 
combination (Burgess 1964) must be included. Currently, the most 
widely used ionisation balance calculations in the literature are 
those of Arnaud & Rothenfiug (1985) and Mazzotta et al. (1998) 
or, for ions of iron, those by Arnaud & Raymond (1992). More 
recently high precision calculations have been made for specific el- 
ements or ions (e.g. Nahar 1995 [Si I]; Nahar 1996 [Si II]; Nahar & 
Pradhan 1997 [C and N]; Nahar 1999 [O]; Nahar, Pradhan & Zhang 
2000 [C, including relativistic effects]). In common with the earlier 
work (Arnaud & Rothenfiug 1985), usually only the zero-density 
limit is considered for recombination rates. 

Burgess (1965) pointed out that non-zero density reduces the 
dielectronic recombination rate because highly excited states ef- 
fectively become part of the continuum. Jordan (1969a) attempted 
to account for this in an approximate way in her ionisation bal- 
ance calculations and Judge et al. (1995) used recombination rates 
corrected for finite density (based on calculations by Summers 
1974a,b) in their analysis of solar data, but none of the more re- 
cent calculations mentioned above address this issue. The effects 
are significant for e Eri, where the transition region pressure (and 
therefore density) is higher than in the Sun - a reduction of the di- 
electronic recombination rate will alter the ionisation balance, af- 
fecting, in particular, the important Li-like and Na-like ions (the 
sensitivity of these ions is apparent from the calculations of Jordan 
1969a). It is well known from solar work that emission measures 
derived from Li-like and Na-like ions often appear anomalous (Bur- 
ton et al. 1971, Dupree 1972, Judge et al. 1995). It has also been 
noted in analyses of other cool stars (e.g. Jordan et al. 1987). Simi- 
lar results are found for e Eri; if the standard Arnaud & Rothenfiug 
(1985) ionisation balance calculations are adopted, the loci of Li- 
like and Na-like ions tend to lie too high compared with those of 
other ions (see Section 4). In view of the sensitivity of these ions 
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to the density-dependent reduction to the recombination rates it is 
necessary to attempt to include this effect in the ionisation balance. 

Summers (1972, 1974a,b) computed recombination rates as 
a function of both temperature and electron density (N e ) for var- 
ious elements, including those of primary interest here (C, N, O, 
Si). Rates computed at electron densities of log /V e [cm ] = 8.0 
and 12 were taken from Summers (1974b). For C, N, O and Si 
these rates were extrapolated linearly against temperature in log- 
log space to find recombination rates appropriate for the mean tran- 
sition region pressure in e Eri (log P e [cm -3 K] = 15.68, Jordan 
et al. 2001a). The recombination rates determined in this way for 
several important ions (C IV, N V and Si IV) are different by sev- 
eral hundred percent from the zero-density limit; this is sufficiently 
large to noticeably affect the ionisation balance. 

These extrapolated recombination rates were used together 
with collisional ionisation rates from Bell et al. (1983) to compute 
the fractional ion populations as a function of temperature for C II 
- IV, N III - V, O III - O VI and Si II - IV. Photoionisation has been 
neglected for these calculations, as has charge exchange for the C, 
N and O ions. Charge exchange (with rates taken from Arnaud & 
Rothenflug 1985) was included for Si. As expected, the ionisation 
fractions computed differ most significantly from those of Arnaud 
& Rothenflug (1985) for Li-like and Na-like ion, but there are also 
differences for several lower ionisation stages. The reduction in re- 
combination rates causes the peak ionisation fraction to occur at 
lower temperatures for many ions (as apparent in Jordan 1969a). 

These effects are illustrated in Fig. 1 which shows ionisation 
balance calculations for C. Two sets of zero-density calculations are 
shown (Arnaud & Rothenflug 1985; Nahar & Pradhan 1997) along 
with Jordan's (1969a) results for a solar pressure and the new com- 
putations for e Eri. Generally, the zero-density calculations (Arnaud 
& Rothenflug 1985 and Nahar & Pradhan 1997) agree well with 
each other but there are significant differences from the finite den- 
sity cases. The agreement between the new calculations and those 
of Jordan (1969a) is good - although Jordan's calculations were in- 
tended for densities in the quiet Sun, she found later that the density 
sensitivity was overestimated and that the results are more applica- 
ble to higher pressures such as those found in active regions, or 
e Eri. The results for the other elements considered (N, O, Si) are 
all similar. The largest differences between the new calculations 
and those by Jordan (1969a) are for Si since the latter did not in- 
clude charge exchange processes. 

The new ionisation balance calculations for C, N, O and Si 
are used in the emission measure modelling (see Section 4) and 
the recombination rates computed for finite density are employed 
in the atomic models of Si and C discussed in Section 5.2. The 
computed ionisation fractions are tabulated in Appendix A. The 
accuracy of these ion fractions is limited by the accuracy of the 
Summers (1972,1974a,b) recombination rates; there is now a real 
need for modern calculations of di-electronic recombination rates 
at finite density. 



4 EMISSION MEASURE DISTRIBUTION (EMD) 
4.1 Emission measure loci 

In this section the fluxes presented in Section 2 are used to con- 
struct loci which constrain the mean emission measure distribution 
Em 3 (T e ). The principles of this method have been discussed at 
length elsewhere (e.g. Jordan & Brown 1981). The notation and 
definitions used by Sim & Jordan (2003b) are adopted here. 



The standard method of constructing emission measure loci 
from line fluxes is only valid for optically thin lines. For optically 
thick lines radiative transfer calculations are required. In order to 
construct a "first-cut" model, however, even optically thick lines 
can be included in the emission measure analysis by assuming that 
they are part of a multiplet which is effectively optically thin. For 
lines in which scattering is the only opacity source this assumption 
should not lead to errors of more than a factor of two. This approach 
is used in this section. In Section 5, full radiative transfer calcula- 
tions are discussed and used to improve the calculations presented 
here. 

Throughout the following discussion, solar photospheric val- 
ues have been adopted for the relative abundances of the elements 
(from Grevesse & Sauval 1998), with absolute values scaled to the 
e Eri metallicity derived by Drake & Smith (1993). The observed 
line fluxes have been converted to fluxes at the stellar surface us- 
ing the stellar parameters adopted by Jordan et al. (2001b, table 1). 
We investigate later the differences which result from adopting the 
recent solar abundances derived by Asplund, Grevesse & Sauval 
(2004). 

4.2 Density-insensitive STIS lines 

For many of the strong emission lines in the STIS wavelength 
range, the contribution function K^e) 1 is almost completely in- 
sensitive to the density. This generally occurs for permitted transi- 
tion in which the upper level is predominantly excited by collisions 
from the ground state. Such lines place the strongest constraints on 
the EMD since their loci are independent of any uncertainty in the 
electron pressure. 

Table 4 lists the density-insensitive lines/multiplets from the 
STIS lines (the insensitivity criterion was taken to be a variation in 
K(T e ) of less than ten per cent resulting from a variation of log P e 
by ±0.3 dex around the mean transition region pressure). Table 4 
also indicates whether the lines may reasonably be assumed to be 
optically thin or whether only the multiplet as a whole is likely to be 
thin (it is stressed that optically thick lines are treated with radiative 
transfer calculations at a later stage, see Section 5). 

The multiplets listed in Table 4 have been used to construct 
emission measure loci. The K(T e ) contribution functions were 
computed with the CHIANTI database (v3.03; Dere at al. 1997; 
Dere et al. 2001). For C, N, O and Si the ionisation fractions were 
computed as described in Section 3. For other elements the Arnaud 
& Rothenflug (1985) ionisation fractions were adopted. The loci 
computed are plotted in Fig. 2. 

The loci suggest a smooth EMD decreasing from the chro- 
mosphere to a minimum somewhere in the mid-transition region at 
around logT e ~ 5.3. There is a sufficiently large number of lines 
to provide temperature coverage throughout this part of the atmo- 
sphere, placing constraints on the EMD that are at least as complete 
as those for the Sun. 

The lightly drawn loci in Fig. 2 (offset by +3 dex for clarity) 
show the constraints derived from the same line fluxes if the Arnaud 
& Rothenflug (1985) ionisation calculations are used for every ion. 
In this case, C IV and Si IV are clearly anomalous by factors > 2 
and N V also appears to lie rather high. However, when the new 
ionisation calculations are adopted (Section 3) the loci from these 

1 An explicit expression for the contribution function (K(T e )) in terms of 
abundance, ionisation fraction, degree of excitation, electron number den- 
sity and Einstein-A coefficient is given by Sim & Jordan (2003b). 
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Figure 1. Ionisation fractions for C (JVj /iVg) as a function of temperature. C I - C V are shown, although only results for C II - C IV are needed in Section 4. 
The curves are labelled with the ionisation stage they represent. The solid curves are the ionisation fractions computed here. The dashed curves are the finite 
density calculations of Jordan (1969a); the dotted curves are the Nahar & Pradhan (1997) ion fractions and the dot-dash curves are from Arnaud & Rothenflug 
(1985). 
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Figure 2. Emission measure loci from density-insensitive lines in the e Eri STIS spectrum. The solid curves are for lines of Si; the dashed curves for C; the 
dotted curve for O; dot-dash curves are used for N, Al and Mg; and the triple-dot-dash curves are for S. The heavy lines show the loci computed using the 
new ionisation fractions from Section 3 and are labelled with the ion they represent. The light lines show the loci computed if all ion fractions are taken from 
Arnaud & Rothenflug. These curves have been off-set by +3 dex, for clarity. The broad, light curve is the mean EMD corresponding to DEMD-A (see Section 
4.5). 
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Table 4. Density-insensitive lines used for emission measure modelling. 
Where an entry is given in parenthesis it is likely that the lines are affected 
by radiative transfer effects which are dealt with in Section 5. 



Ion 


A (A) 


Assumed opacity 


(C II 


1335 


multiplet effectively thin) 


C III 


1175 


thin 


C IV 


1550 


thin 


N V 


1240 


thin 


v 


1371 


thin 


(Mgll 


2800 


multiplet effectively thin) 


Aim 


1854 


effectively thin 


Si II 


1810 


multiplet effectively thin 


Si II 


1530 


multiplet effectively thin 


Si II 


1308 


multiplet effectively thin 


Si II 


1260 


multiplet effectively thin 


Si II 


1195 


multiplet effectively thin 


Si III 


1206 


effectively thin 


Si IV 


1400 


thin 


SO 


1255 


multiplet effectively thin 


S IV 


1406 


thin 



Li-like and Na-like ions are in much better agreement with those of 
other ions. The agreement is still not perfect, in particular the Si IV 
locus still lies higher than that of C III by about 30 per cent, but this 
difference is small compared to the likely residual uncertainty in the 
ionisation balance calculations. The Al III locus is also rather high; 
although this may be the result of blending in the observed Al III 
line, it may also be the result of errors in the ionisation balance; 
like C IV, Al III is an alkali but no new ionisation calculations have 
been performed for this ion. 

It is concluded that all the loci which appear significantly 
anomalous when the Arnaud & Rothenflug (1985) ionisation bal- 
ance calculations are adopted are likely to be the result of errors 
in the recombination rates used in that ionisation balance, owing 
to the density dependence of dielectronic recombination. This is in 
contrast to the analysis of solar data by Judge et al. (1995); they 
find that loci derived for alkali-like ions are still discordant, even 
if they account for the density-sensitivity of dielectronic recombi- 
nation and they interpret this as evidence of dynamic or diffusive 
effects (although they do state that uncertainties in the recombi- 
nation rates may still be significant). 2 This issue cannot be fully 
resolved without reliable recombination rates. 

4.3 Density-sensitive STIS lines 

The density-insensitive STIS loci can be supplemented by loci 
derived from density-sensitive lines/multiplets. Table 5 lists the 
density-sensitive multiplets from the STIS data. Most of these mul- 
tiplets and the associated atomic data have been discussed by Jor- 

2 However, Judge et al. (1995) used a mixture of lines above and below the 
H Lyman edge at 911 A, whereas all the lines used here lie above 911 A. 
Judge et al. (1995) rejected absorption by the Lyman continuum as a source 
of the discrepancies that they found, partly on the basis that centre-to-limb 
differences have not been observed. However, Burton et al. (1973) estab- 
lished the magnitude of this by observing the limb-to-disc ratio of pairs of 
optically thin lines from a given ion, where one was found at A > 911 A and 
the other at A < 911 A. This amounts to a factor of up to ~ 3 between 
A ~ 765 and 835 A. These ratios are independent of any atomic physics, 
abundances or intensity calibration. The effect expected for stellar observa- 
tions should, of course, be less than this. 



Table 5. Density-sensitive lines used for emission measure modelling. 



Ion 


A (A) 


Opacity 


NIV 


1486 


thin 


O III 


1663 


thin 


O IV 


1400 


thin 


O V 


1218 


thin 


Si III 


1892 


thin 


Sim 


1300 


thin 



dan et al. (2001a) where they were used to deduce the mean tran- 
sition region pressure. The Si II 2334-A multiplet is not used here; 
this multiplet is density-sensitive but there are line blends and ob- 
vious problems with the atomic data. 

The loci deduced from the density-sensitive loci are plotted 
in Fig 3. In general there is good agreement between the density- 
sensitive and density-insensitive loci when logP e [cm~ 3 K] = 
15.68 is adopted. The N IV and, in particular, the O IV loci lie 
lower than those from density-insensitive lines at similar tempera- 
tures. A similar effect has been noted in the Sun; Judge & Brekke 
(1994) found that loci derived from density-sensitive lines lie lower 
than those from resonance lines by a significant factor. In part this 
may be due to uncertainties in the atomic data for these transitions 
(it is known that there are inadequacies in the O IV atomic data - 
see the discussion by Jordan et al. 2001a). There are also compli- 
cations in the interpretation of the loci from density-sensitive mul- 
tiplets resulting from inhomogeneity in the atmosphere (Judge & 
Brekke 1994). If the atmosphere is inhomogeneous then loci de- 
rived from density-insensitive lines will represent unbiased aver- 
ages over the inhomogeneous structure. In contrast, the density- 
sensitive lines will form predominantly in regions with higher or 
lower density than average resulting in a pressure-weighted aver- 
age. To investigate this effect it is necessary to quantify the nature 
of the sensitivity of the density-sensitive lines. Fig. 4 shows loci de- 
rived for the density-sensitive lines at different pressures. For most 
of the lines, a higher pressure means that less flux is formed for a 
given emission measure (i.e. the emission measure locus lies higher 
for a given flux) and vice versa. The Si III 1300- A multiplet is an 
exception because of the influence of the density on the metastable 
3 P populations. The N IV, O III and Si III 1892-A loci are the most 
sensitive to errors in the pressure, changing by about 50 per cent 
when the pressure is varied by a factor of two. O IV and O V are 
interesting because, although less sensitive than O III, they respond 
more to a high density than a low density. Thus if there is a range 
of densities present in the atmosphere those places with density 
lower than average will not contribute significantly more than aver- 
age, but those regions with high density will contribute significantly 
less. This may partially explain why the O IV locus appears to lie 
too low. 

4.4 FUSE lines 

The multiplets measured in the FUSE spectra provide useful sup- 
plementary constraints. Table 6 lists the multiplets observed with 
FUSE that have been used to construct emission measure loci. The 
loci are shown in Fig. 5. 

Most of the FUSE loci agree well with each other, and with 
the STIS loci. The FUSE S IV locus is not consistent with the STIS 
loci. The S IV measurement error is large (~ 40 per cent) but not 
sufficient to explain the discrepancy. The most probably explana- 
tion is a systematic error in the S IV atomic data. 
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Figure 3. Emission measure loci from density-sensitive lines in the e Eri STIS spectrum. The solid curves are for lines of Si; the dotted curve for O; the 
dot-dash curve for N. The light broad curve shows the same EMD plotted in Fig. 2. 
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Figure 4. Emission measure loci from density-sensitive lines in the e Eri STIS spectrum. The curves are for the same ions as indicated in Fig 3. The solid 
curves were calculated using the best mean electron pressure, the dotted curves using a pressure smaller than the mean by a factor of two and the dashed curves 
using a pressure higher than the mean by a factor of two. Note that the scale is significantly expanded compared to Figs 2 and 3. 



4.5 Preliminary mean EMD 

The procedure of constructing emission measure distributions from 
emission measure loci has come under close scrutiny in recent years 
(e.g. Judge, Hubeny & Brown 1997; Mcintosh, Brown & Judge 
1998). Since the EMD is a continuous function but we have only 
a finite number of constraints, each with an associated uncertainty, 
it is not possible to determine a unique EMD from the data alone. 



Thus the EMD obtained will depend on how it is obtained. How- 
ever, as stressed by Judge et al. (1997), the EMD determined to be 
consistent with the data is still important since it represents the most 
information that can be extracted from the data alone, and the range 
of distributions that are consistent with the data is an indication of 
the completeness of the information provided by the data. 

Here, the mean EMD (Em 03 (T e )) has been constructed 
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Figure 5. Emission measure loci from lines in the e Eri FUSE spectrum. The solid curve is for the Si multiplet; the dashed curves for C; the dotted curve for 
O; the dot-dash curve for N; the triplet-dot-dash curves for S. The light, broad curve shows the mean EMD from Fig. 2 for comparison. 



Table 6. FUSE lines used for emission measure modelling. The last column 
gives the FUSE channel in which the lines are observed. Only the pointing 
of the 1-LiF is accurate. The other channels may lose photons due to channel 
drift making the measurements from these channels prone to large errors. 



Ion 


A (A) 


Opacity 


Channel 


Cll 


1010 


effectively thin 


lA-LiF 


CII 


1037 


(effectively thin) 


lA-LiF 


C III 


1175 


thin 


lB-LiF 


C III 


977 


effectively thin 


lB-SiC 


NlII 


991 


thin 


lA-LiF 


O VI 


1035 


thin 


lA-LiF 


Si III 


1110 


thin 


lB-LiF 


S IV 


1070 


thin 


lA-LiF 


S VI 


933 


thin 


lB-SiC 



by first obtaining the differential emission measure distribution 
(DEMD) denned by 

(see Sim & Jordan 2003b for definitions of the quantities in this and 
the following equation). A first-cut DEMD was used to compute 
fluxes for each of observed multiplets using the integral 

K = / TC) dTof?; dlogre • (2) 

These were compared with the observed fluxes. The DEMD was 
then modified iteratively until reasonable agreement (~ 50 per 
cent) was obtained between the fluxes predicted by the DEMD and 
the observed fluxes. Greatest weighting was given to the density- 
insensitive lines. During this process, it was assumed that the 
DEMD remains smooth on scales A log T e [K] < 0.1 dex, which is 



comparable to the typical scale of overlap between adjacent emis- 
sion measure loci from the STIS data. 

The DEMD was only specified over the temperature range 
3.8 < logT e [K] < 5.5 since outside this range there are insuffi- 
cient lines to provide robust constraints. This means that lines that 
form partially outside this temperature range (in particular those of 
O VI and O V) will be underpredicted by the DEMD. The resulting 
DEMD (DEMD- A) was used to compute a mean EMD (EMD- A) 
which is overplotted in Figs. 2, 3 and 5. The fluxes computed from 
DEMD-A using equation (2) are compared with observed fluxes in 
Table 7. 

DEMD-A predicts fluxes for most of the high-temperature, 
optically thin lines (C III, C IV, N V, Si III, Si IV, S IV, N IV, 
O III) that are in agreement with the observations to within ~ 30 
per cent, which is comparable to the likely measurement errors for 
the weak lines, and the uncertainty in the ionisation balance for 
the alkali-like ions. Among the high-temperature lines, the poorest 
agreement is for O IV 1400 A and S IV 1070 A. Both of these mul- 
tiplets have already been identified as anomalous from the EMD 
loci discussed above. Note that for both Al III and S VI the under- 
prediction is probably the result of errors in the ionisation balance 
- these are both alkali-like ions for which the finite density has not 
been accounted for in the recombination rates (see Section 3). 

Among the lower-temperature ions (C II, Mg II and Si II) the 
agreement is less good. However, these multiplets are likely to be 
affected by opacity (making equation 2 invalid) - their discussion 
is postponed to Section 5 where full radiative transfer calculations 
are presented. 



5 THE PRELIMINARY ATMOSPHERIC MODEL 
5.1 Construction of the model 

In this and subsequent sections, EMDs are used to construct new 
atmospheric models for e Eri. These models are used to perform ra- 
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Table 7. Comparison of calculated fluxes (F c ) and observed fluxes (F a ) 
for emission measure distribution DEMD-A. The error column gives the 
percentage error on the measured flux. The top part of the table gives the 
density-insensitive STIS lines, the middle part gives the density-sensitive 
STIS lines and the lower part gives the FUSE lines. 



Ion 


A (A) 


% error 


(F c - F )/F x 100% 


CII 


1335 


3 


-7 


C III 


1175 


10 


13 


C iv 


1550 


2 


-19 


N V 


1240 


4 


-29 


O V 


1371 


5 


-40 




9 8 no 


2 


-41 


Aim 


1854 


12 


-60 


Sill 


1810 


1 


-55 


Si ii 


i s^o 

1JJU 


i 


-OO 


Si II 


1308 


5 


-80 


Si II 


1 960 


1 9 




Si II 


1195 


30 


-71 


Si III 


1206 


3 


-17 


Si IV 


1400 


3 


-29 


S II 


1255 


10 


_4 


*\ TV 




90 


Q 


NIV 


1486 


33 


30 


O III 


1663 


17 


16 


O IV 


1400 


6 


89 


O v 


1218 


6 


-22 


Sim 


1892 


4 


5 


Sim 


1300 


7 


-22 


Co 


1010 


25 


30 


Cn 


1037 


6 


56 


C m 


977 


2 


52 


C m 


1175 


2 


11 


Nm 


991 


24 


35 


O VI 


1035 


4 


-50 


Sim 


1110 


30 


12 


S IV 


1070 


43 


125 


S vi 


933 


33 


-68 



diative transfer calculations for the important optically thick lines, 
the results of which are then discussed. 

The atmospheric models presented are all constructed in one- 
dimension under the assumption of hydrostatic equilibrium, in- 
cluding both thermal and microturbulent contributions to the pres- 
sure. Full details of the method of constructing atmospheric mod- 
els from EMDs are presented elsewhere (see e.g. Griffiths & Jor- 
dan 1998 and Harper 1992 for all the relevant background or Sim 
2002 for a complete discussion of the iteration scheme used here). 
The Thatcher et al. (1991) model is used to describe the atmo- 
sphere below the temperature range to which the DEMD extends 
(logT e [K] < 3.8). 

An important ingredient in the modelling is the microturbu- 
lent velocity in the chromosphere/transition region. Sim & Jordan 
(2003a) used line width measurements from the e Eri ultraviolet 
data to determine the most-probable non-thermal speed (£(T e )) in 
the upper atmosphere. For the models presented here, these results 
have been fitted with the simple form 

c/r\-f -31.56 + 10.716 log T e 3.8 < logT e < 4.9 \ 
^ U >-\ 20.95 logT e > 4.9 J 

(3) 



which is an adequate representation of the line-width measure- 
ments. Initially, it is assumed that all the observed non-thermal mo- 
tion can be classified as microturbulence (i.e. £(T e ) m i C ro = £(T e )). 
Later, a small amount of macroturbulence is introduced (Section 6). 
For logT e [K] < 3.8, the microturbulence is extrapolated linearly 
in log T e to agree with the measured photospheric microturbulent 
velocity (1.25 km s _1 ; Drake & Smith 1993) at the temperature 
minimum. 



5.2 Radiative transfer calculations 

Radiative transfer calculations were performed with the MULTI 
code (Scharmer & Carlsson 1985a,b; Carlsson 1986). In addition 
to the atmospheric models, as input MULTI requires atomic mod- 
els for the ions that give rise to the lines under consideration. Brief 
details of the adopted atomic models are given below. 



5.2.1 Hydrogen 

The H model used contains levels n = 1 to n = 8 and the con- 
tinuum. The radiative transfer calculations for H make use of the 
Hubeny & Lites (1995) modifications to MULTI which allow the 
effects of Partial Redistribution (PRD) to be taken into account for 
Lyman a and (3. (See Hubeny & Lites 1995 and Sim 2001 for dis- 
cussions of PRD effects in cool star atmospheres.) 



5.2.2 Magnesium 

A ten level Mg model, developed by P. Judge, was used. The model 
describes the first three ionisation stages of Mg. PRD is included 
in the h and k-lines of Mg II using the Uitenbroek (1989) modifica- 
tions to MULTI. 



5.2.3 Silicon 

In view of the large number of Si lines observed in the STIS data, a 
new 36 level model of Si has been constructed. This model focuses 
on Si II - IV and contains the ground states of Si I and Si V, for the 
sake of the ionisation calculations. The metastable 3p 2 1 D, X S and 
3s3p 3 5 S levels of Si I are also included. 51 radiative transitions 
are included in the model (all treated with Complete Redistribu- 
tion [CRD] only). The atomic data for these transitions (gf-values 
and rates for collisional excitation by free electrons) were extracted 
from the CHIANTI database. Recombination rates (radiative plus 
dielectronic) were taken from Summers (1974a,b, see Section 2) 
and photoionisation rates were obtained from the on-line Opacity 
Package database (Cunto et al. 1993). Collisional ionisation rates 
and charge exchange rates for collisions with both hydrogen and 
helium were taken from Arnaud & Rothenflug (1985). 



5.2.4 Carbon 

A new 33 level model for C has been constructed to complement 
the new Si model. Like the Si model it is focused on the second to 
fourth ionisation stages, containing only the ground states of the 
first and fifth stages for the ionisation calculations. 52 radiative 
transitions are included in the model (all treated with CRD). The 
atomic data were drawn from the same sources as used for the Si 
model. 
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Table 8. Fluxes computed using MULTI (F c ) compared with observed fluxes (F ) for lines of H, Mg, Si and C. F c is given in units of 10 -13 ergs cm~ 2 s _1 . 
The first section gives the density- insensitive lines from the STIS spectrum, the second gives the density-sensitive lines from STIS and the third gives the lines 
from FUSE. Computed fluxes are given for models A, B, C and D. Following the model B fluxes, values are given for computations from a model identical to 
model B except that the Asplund et al. (2004) relative abundances are adopted. 



Ion 


Ao 
(A) 


Model A 


F c 

Model B 


Model C 


Model D 


Model A 


(F c - F )/F 
Model B 


x 100% 
Model C 


Model D 


rl 1 


m c ao i 
1215.6 / 1 


A1A 
034. 


A A 1 f A\(W 
441 (410) 


302 


OOA. 
220 


1 O/l 

124 


1 OA f1 CA) 
1 /0 (156) 


88.8 


io C 

3 / .5 


C II 


1 i i a en 
1 334.532 


1 in 

3.30 


i a 1 /o /in\ 
3.41 (2.49) 


1 ai 
1.63 


1 oc 
3.25 


oc c 
/5.5 


01 A /"30 A \ 

81.4 (32.4) 


111 
-13.3 


oo n 
/2.9 


L. II 


1 11C AO Ul 

1335.68 bl. 


8.zy 


o Ac /A oa\ 
8.68 (6.26) 


A f\A 

4.04 


8.81 


O 1 1 

/1. 3 


on i /on i\ 
/9.3 (29.3) 


1 A C 

-16.5 


oo n 
82.0 


C III 


1 1 O/l C\11 


0.341 


n cc/: /n /in/i\ 
0.556 (0.404) 


0.453 


n con 
0.589 


10 n 
32.0 


1 C A / 1 C\ 

35.6 (-1.5) 


1 n c 
10.5 


A 1 O 

43. / 


C III 


I 1 OC OA1 

I I /5.263 


0.423 


n /in /n a 1 q\ 
0.43 / (0.318) 


n ico 
0.35 / 


n a A 1 
0.44 1 


1 "7 

3. / 


A A / OO A \ 

6.6 (-22.4) 


1 o n 
-12.9 


O A 

/.6 


C in 


1 nc An 1,1 

1175.69 bl. 


1.908 


i n£ i/1/i on\ 

1.961 (1.429) 


1.602 


1.976 


28.1 


31.6 (-4.1) 


7.5 


32.6 


C in 


1 175.987 


0.421 


A /lOO /A Q 1 C\ 

0.433 (0.315) 


0.353 


0.438 


2.7 


C A / 01 O) 

5.6 (-23.2) 


-13.9 


6.8 


C III 


I 1 oa ion 

I I /O.J /O 


0.32o 


n c /i i /n incA 
0.543 (0.395) 


n A A A 
0.444 


a ck 

0.545 


14.8 


ion / 1y1 1\ 

18.0 (-14.1 ) 


-3.5 


1 c 
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1548.187 


3.83 


A in /i 1 a\ 

4.39 (3.19) 


4.32 


4.36 


-31.5 


O 1 C / /1 A\ 

-21.5 (-42.9) 


-22.7 


-22.0 


C IV 


1 CCA 7T1 

1550.772 


1.92 


o o a / 1 c n \ 

2.20 (1.59) 


2.16 


2.18 


-31.9 


OOA/ /lO A \ 

-22.0 (-43.6) 


-23.4 
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Mg II 


OOnA "J CO 
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Q /I O /I 1 0\ 
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3.0 
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15.1 


-30.2 
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2803.531 


i n< 
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11 O / 1 (11) 
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ooo 
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1 c c 
-15.5 


An/ 1 A O) 

-6.0 (-16.8) 


1 o 
-Y.I 


/in c 
-40.5 


Si II 


1 1 nn A 1 A 

i iyu.4io 


n m i 
0.031 


n nAn /n accn 
0.060 (0.055) 


n r\~i A 
0.0/4 


n ncA 
0.056 


An n 
-69.0 


/inn/ /i c n\ 
-40.0 (-45.0) 


oa n 
-26.0 


a a n 
-44.0 


Si II 


1193.290 


0.033 


0.057 (0.054) 


0.075 


0.043 


-63.3 


-36.7 (-40.0) 


-16.7 


-52.2 


Si II 


1194.500 


0.092 


0.131 (0.122) 


0.170 


0.101 


-48.9 


-27.2 (-32.2) 


-5.6 


-43.9 


Si II 


1 1 no an/i 
I IV /.jV4 
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0.013 


n na/i /n mi) 
0.034 (0.033) 


n t\A c 
0.045 


n noo 
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A TIC /A 1 ATI 

0.215 (0.192) 
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-19.4 


12.6 (0.5) 
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Si II 


1 OA/1 "71 O 
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0.305 
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Si II 


n/c nr\i 
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0.107 
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Si II 
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-43.0 (-30.4) 
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-46.6 


Si II 


1309.276 


0.071 


0.107 (0.095) 


0.126 


0.096 


-68.0 


-51.8 (-57.2) 


-43.2 


-56.8 


Si II 


1526.707 


0.294 


0.570 (0.508) 


0.661 


0.478 


-42.1 


12.2 (0.0) 


30.1 


-5.9 


Si II 


1 533.431 
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4.03 
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37.0 (0.0) 


-33.3 


44.4 


CII 


1010.371 


0.055 


0.055 (0.040) 
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0.453 


0.589 


24.4 


27.8 (-7.1) 


4.1 


35.4 


C III 


1175.263 


0.425 


0.437 (0.318) 


0.357 


0.441 


0.0 


2.8 (-25.2) 


-16.0 


3.8 


C III 


1175.69 bl. 


1.908 


1.961 (1.429) 


1.602 


1.976 


25.6 


29.1 (-5.9) 


5.5 


30.1 


C III 


1175.987 


0.421 


0.433 (0.315) 
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a Line affected by interstellar absorption. 
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5.3 Model A 

Using the assumptions mentioned in Section 5.1, DEMD-A was 
used to compute a new atmospheric model for the chromo- 
sphere/lower transition region. This model was used by Sim & Jor- 
dan (2003a); the most important atmospheric parameters (T e , P e , 
Ph) are shown in their fig. 3. 

Together with the atomic models discussed in Section 5.2, 
Model A was used to compute fluxes for the observed lines of H I, 
C II - IV, Mg II and Si II - IV. These fluxes are given in Table 8, 
where they are compared with the observed fluxes. For comparison, 
the run of electron temperature against mass column density from 
the Thatcher et al. (1991) model was also used to compute fluxes - 
these are given in Table 9. Since the Thatcher et al. (1991) model 
was developed using chromospheric lines, only lines of relatively 
low temperature ions are listed in Table 9 - for all of the higher- 
temperature lines (which were outside the scope of their study) the 
Thatcher et al. (1991) model substantially overpredicts the fluxes. 

In general, model A can reproduce the observed fluxes of 
Mg II, C III - IV and Si III - IV within an accuracy of several tens 
per cent (see Table 8). The poorest agreement for lines of these ions 
is for C III 977 A, the observed flux for which may be subject to 
significant uncertainty due to the correction for channel drift in the 
FUSE spectrum. 

The greatest failings of model A relate to C II and Si II; 
although the weak C II 1010- A multiplet is predicted to within 
observational errors (almost), the strong resonance multiplets 
(1335 A and 1037 A) are overpredicted by factors of ~ 2. These 
multiplets are predicted to be significantly stronger by the radiative 
transfer calculations (Table 8) than the optically thin approxima- 
tion (Table 7) - this is due to scattering in the lines (see Section 
5.3.3). Although interstellar absorption affects C II 1334.5 A and 
1036.3 A, it cannot be responsible for the discrepancies in the 
stronger 1335.7-A and 1037.0-A lines. In contrast, the strong Si II 
lines are all underpredicted by as much as 70 per cent. Although 
the Thatcher et al. (1991) model does slightly better than model A 
for some of the chromospheric Si II lines, it does much less well 
for the hotter transition region lines (i.e. C II) - as noted above, 
the Thatcher et al. (1991) model was not developed to address the 
transition region in detail. 

While small alterations to the DEMD could improve the 
agreement amongst the higher temperature lines, it is not possi- 
ble to simply attribute these results to failings in the determina- 
tion of the DEMD. By examining the emission measure loci in 
Figs. 2 and 5, it is clear that the C II lines form over a wide range 
4.0 < logT e (K) < 4.5. To reduce the predicted C II resonance 
line fluxes would require the DEMD to be reduced in this tem- 
perature range. However, this cannot be done without reducing the 
predicted fluxes of Si III or Si II which form in the same tem- 
perature range. Such a reduction is unacceptable since the Si III 
fluxes are not overpredicted by model A, and the Si II fluxes are 
already underpredicted to a significant extent. That simple modifi- 
cations to the DEMD cannot account for the C Il/Si II discrepancy 
was checked by constructing ~ 20 different models with modified 
DEMDs, none of which were able to significantly improve the pre- 
dicted C II fluxes without substantial adverse effects on the Si II 
and/or III fluxes. 

Therefore, an alternative explanation is required for the 
C Il/Si II discrepancy. Potential explanations are discussed below 
and used in Section 5.4 to construct two more models for the outer 
atmosphere of e Eri (models B and C). 



Table 9. Fluxes computed using MULTI (F c ) and the Thatcher et al. (1991) 
model compared with observed fluxes (F ) for lines of H, Mg, Si II and 
C II. F c is given in units of 10 -13 ergs cm -2 s _1 . The upper section gives 
the lines from the STIS spectrum and the lower gives the those from the 
FUSE spectrum. 



Ion 


Ao 

(A) 


F c 


(F c - F )/F x 100% 


H l a 


1215 671 


1 175 


639 


Cll a 


1334.532 


67.9 


3512 


CII 


1335.68 bl. 


213.9 
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Mg II a 


97Q6 ^S9 


jyo 


jU.j 


Mg II a 


ZoUj.JjI 




1 n ^ 


Si II 


1 190 416 
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<si TT 
ol 11 


I 1 9QO 

I I y3.Z.yU 


U.JO 1 


^9^ 


Si II 


1 1 Q4 son 


0.829 
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<si TT 
Ol 11 


1 1 Q7 


O 9S6 


1 ^fi 

IJD 


Sill a 


1260.422 


1.02 


434 


Sill 


1264.738 


2.03 


302 


Sill 


1265.002 


0.470 


173 


Sill 


1304.370 


0.303 


128 


Sill 


1309.276 


0.455 


105 


Sill 


1526.707 


1.07 


111 


Si II 


1533.431 


1.53 


152 


Sill 


1808.013 


6.26 


-2.0 


Si II 


1816.928 


11.75 


12.8 


Sill 


1817.451 


2.89 


-9.7 


Hi" 


972.537 


9.70 


1440 


Hi" 


1025.722 


29.0 


692 


CII 


1009.858 


1.32 


10900 


CII 


1010.083 


2.25 


8233 


CII 


1010.371 


2.98 


6378 


Cll a 


1036.337 


19.4 


16067 


CII 


1037.018 


22.1 


7495 



a Line affected by interstellar absorption. 



5.3.1 Atomic data for Si 11 

The limitations of the atomic data used in the radiative transfer cal- 
culations are almost certainly partially responsible for the apparent 
discrepancy amongst lines of Si II. 

The most important uncertainties are in the electron collision 
rates. The rates used in the Si II model are from Dufton & Kingston 
(1991). Although it is difficult to asses the accuracy of collision 
rates, Dufton & Kingston (1991) suggest that their rates should be 
accurate to at least ±20 per cent. Their calculations give rates be- 
tween the ground 2 P levels and all the excited states used in the Si II 
model, and between the metastable 4 P and excited 3s3p 2 2 D levels 
(which give rise to the 1810-A multiplet), but no collision rates are 
provided from the 4 P levels to other excited states. The absence of 
excitation rates from the metastable 4 P states may mean that the 
populations of all the excited states above 3s3p 2 2 D are underesti- 
mated - this would cause the 1530, 1260, 1309 and 1190-A mul- 
tiplets to be systematically underpredicted. Assuming that the col- 
lision strengths from the 4 P states to levels higher than 3s3p 2 2 D 
are the same as those to the 3s3p 2 2 D 3 / 2 level, it is found that the 
predicted fluxes in the multiplets mentioned above are increased 
by 53, 25, 73 and 32 per cent, respectively. Although only a crude 
estimate, this suggests that excitation from 4 P could be important 
and may explain why many of the Si II transitions are significantly 
underestimated by the models. 

In addition to direct excitation from the ground state, the 4s 2 S 
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level population is significantly affected by collisional excitation to 
4p 2 P followed by radiative decay (Jordan 1969b). Thus the com- 
puted 1530- A flux is probably subject to the largest uncertainties 
due to the atomic data. The 1810-A flux is subject to the smallest 
uncertainties, being dominated by the errors in the collision rates 
between the low-lying levels (at most a few tens per cent). (See 
also Jordan et al. 2001a.) 



5.3.2 Atomic data for C II 

The collision rates for C II are from Blum & Pradhan (1992) who 
give rates between all levels that are relevant to the formation of the 
1335 and 1037-A multiplets. The authors suggested that these rates 
are accurate to between 20 per cent and a factor of 2, with greater 
uncertainty for higher levels. This suggests that the data should be 
better for the 1335-A multiplet than for the 1037-A multiplet. 

5.3.3 Geometric complications 

The C II resonance line fluxes are computed to be significantly 
larger when radiative transfer is taken into account, compared to 
the optically thin limit (compare Tables 7 and 8). This is because 
the calculations in Section 4.5 and Table 7 assume that the frac- 
tion of photons emitted in a particular line which escape the atmo- 
sphere is ~ 1/2, as given by the usual geometric dilution factor 
(the optically thin limit). For scattering dominated optically thick 
lines (such as the C II resonance lines) that assumption underes- 
timates the surface flux since scattering turns around many of the 
photons that are initially emitted downwards into the atmosphere 
so that they eventually emerge through the upper boundary. For a 
plane-parallel layer this backscattering effect can increase the flux 
by up to a factor of two (for a comprehensive discussion of this 
and related effects in the solar case, see Pietarila & Judge 2004). In 
the MULTI calculations presented here, the intensity is computed 
at the surface of the star as a function of both wavelength and angle 
accounting for radiative transfer effects (including backscattering) 
and this is integrated over angles to compute the flux. Thus for an 
homogeneous, plane-parallel atmosphere, the results presented in 
Section 5.2 and Table 8 are more reliable than those given in Sec- 
tion 4.5. 

In a real stellar atmosphere, however, the situation is more 
complex; assuming that t Eri is like the Sun, the transition region 
emission will mostly originate in the bright supergranulation cell 
boundaries which occupy only a fraction of the stellar surface (see, 
for example the solar C II lines discussed and modelled by Judge, 
Carlsson & Stein 2003). We have derived a filling factor of emit- 
ting material in the upper transition region of e Eri of ~ 0.2 (Sim 
& Jordan 2003b). In such a geometry, much of the emission will 
not travel solely within the cell boundary material but will escape 
sideways into the cell interiors which are likely to be hotter and 
less dense. This will alter the probability that the emission is sub- 
sequently scattered upward to escape from the star. Thus a single- 
component model may be expected to predict erroneously large 
fluxes for lines in which there is significant scattering opacity. (See 
also Jordan, Smith & Houdebine 2005.) 

The simplest test for this interpretation of the overprediction 
of the C II resonance line fluxes is to examine C II lines that are 
not affected by back-scattering. Unfortunately the ideal lines for 
this examination (the optically thin C II 2s 2 2p 2 P - 2s2p 2 4 P lines) 
are too weak to measure in the STIS spectrum and all the model 
predictions are consistent with this non-detection. The only lines 



that can be used for this test are the weak C II 1010-A lines in 
the FUSE spectrum; the strongest of these lines is measured to an 
accuracy of 15 per cent. Since they originate in a relatively high- 
lying energy level, the atomic data (collision rates) for these lines 
are likely to be significantly poorer than for the other C II lines dis- 
cussed here. Nevertheless, these lines are not significantly overpre- 
dicted by model A (the strongest line agrees with the observations 
to within 20 per cent). Since the lines of the 1010-A multiplet are 
not affected by back-scattering this supports the argument that the 
discrepancy between the observed and computed C II resonance 
line fluxes originates in the treatment of scattering in such lines. 

In order to test the hypothesis that the overprediction of the 
C II resonance lines is due to the simplistic 1-D geometry, multi- 
component models need to be constructed and used to perform 
multi-dimensional radiative transfer calculations. Such work lies 
beyond the scope of this paper. In Section 6 a complete atmo- 
spheric model (model B) is constructed under the hypothesis that 
the C Il/Si II discrepancy is due to errors in the computed C II 
resonance line fluxes. This model has been used to perform a test 
(based on the observed CI I line) which probes the C II resonance 
line strength within the atmosphere. That test supports the evidence 
from the 1010-A multiplet that the discrepancy lies in the compu- 
tation of the C II resonance line surface fluxes. 

5.3.4 Abundances 

An alternative explanation for the C Il/Si II discrepancy is an error 
in the assumed abundances. As mentioned in Section 4.1, the rel- 
ative abundances adopted here are based on the solar values given 
by Grevesse & Sauval (1998). Recently, Asplund et al. (2004) pre- 
sented a new set of solar abundances which differ in the logarithm 
from those of Grevesse & Sauval (1998) by -0.13, -0.14, -0.17 and 
-0.04, for C/H, N/H, O/H and Si/H, respectively. Given that the rel- 
ative changes of C, N and O are very similar (factors of 1.35 to 
1.48) they will have little influence on the shape of the EMD de- 
duced based on lines of these elements - if these new abundances 
were adopted, a small global increase in the EMD (equivalent to 
a small reduction in the adopted metallicity of e Eri) would lead 
to computed fluxes which differ only very slightly from those ob- 
tained using the earlier abundances. Adopting the new abundances, 
however, would alter the relative abundance of C to Si by a factor 
of 1.23. Although this factor would act in the correct sense to ad- 
dress the C Il/Si II problem here, it is insufficient to explain most of 
the discrepancy. Therefore, the standard Grevesse & Sauval (1998) 
abundances will be retained throughout this investigation (although 
the influence of adopting the Asplund et al. 2004 abundances will 
be quantified in Section 6.3). 

To bring the C II resonance lines and Si lines into agreement, 
either the C abundance could be made smaller or the Si abundance 
made larger. If the assumed abundances were wrong, one would 
expect to see a discrepancy between the lines of higher ions (C III, 
IV or Si IV) and lines of other species. Unfortunately, neither C IV 
nor Si IV are ideal for comparison since they are very sensitive to 
the ionisation balance calculation (Section 3). The observed O III 
lines form at similar temperatures to the C III 1 175-A multiplet and 
the Si IV lines. Based on the multiplet fluxes in Table 7, there is 
no obvious discrepancy between O III and C III. Si IV is differ- 
ent, however, being underpredicted rather than overpredicted; thus 
it is concluded that if abundance variations are responsible for the 
C Il/Si II it is the Si abundance that is likely to be in error. 

In order to remove the C Il/Si II discrepancy, the Si abun- 
dance needs to be increased by about a factor of 2.4, relative to 
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the value adopted in model A. The interpretation of such an in- 
crease requires some care. Abia et al. (1988) have measured the 
photospheric Si abundance in e Eri. They found the logarithmic Si 
abundance relative to the Sun to be [Si/H] = —0.11 dex. Based 
on their discussion of likely errors, the uncertainty in this value 
should not be more than ±0.25 dex. More recent studies have 
found similar results: [Si/H] = —0.10 and —0.04 (Bodaghee et 
al. 2003; Zhao et al. 2002, respectively) with much smaller errors 
(~ 0.05 dex). These measurements are all consistent with the as- 
sumed metal abundances in model A which are based on Drake & 
Smith (1993). It suggests, therefore, that a proposed increase in the 
transition region Si abundance would be an increase relative to the 
stellar photospheric abundance. Note that Zhao et al. (2002) de- 
duce a C abundance which is slightly higher than solar and thus, if 
adopted, would make the discrepancy worse. However, since their 
C analysis involves only one spectral line the error in this result 
may be significantly larger than that in the Si abundance. 

It is concluded that if abundance variations are responsible for 
the C Il/Si II discrepancy, the Si abundance enhancement required 
is relative to the e Eri photospheric abundance. There is consider- 
able interest in the relative abundance of elements in the transition 
regions and coronae of cool stars. Many authors (see for exam- 
ple Jordan et al. 1998 and references therein), have discussed the 
possibility that the relative abundances of elements in the corona 
are different from those in the photosphere to a degree which cor- 
relates with the First lonisation Potential (FIP) of the element. 
It is suggested that, in the upper atmosphere, elements with high 
FIPs (> 10 eV) have lower abundances relative to those with low 
FIPs. The first evidence of this "FlP-effect" was found in the so- 
lar wind (Crawford, Price & Sullivan 1972). The physical origin of 
the effect is unknown. Based on studies of the solar atmosphere, 
Feldman et al. (1992) give coronal elemental abundances, assum- 
ing that the FlP-effect acts to increase the low-FIP element abun- 
dances. They give a relative coronal abundance for Si to C which is 
a factor 3.8 greater than the photospheric relative abundance. Lam- 
ing, Drake & Widing (1995) used full disk solar observations to 
confirm the coronal abundances given by Feldman et al. (1992). 
Their work found evidence of the FlP-effect in the corona, but they 
found no suggestion of a FlP-effect operating at lower temperatures 
(log T e (K) < 6). Jordan, Macpherson & Smith (2001) studied spa- 
tial variations in the Si III/C II intensity ratio across regions of the 
quiet Sun and a coronal hole. They found no differences in excess 
of 0.2 dex. Since the FIP effect is not thought to operate in coro- 
nal holes this suggests that the FIP enhancement of Si in the quiet 
Sun is less than a factor of 1.6. These results are surprising. Since 
the first ionisation for most elements takes place at chromospheric 
temperatures, one would expect an effect associated with the first 
ionisation potential to originate in the chromosphere and thus to be 
apparent in the transition region, but the solar work suggests that 
the effect only occurs at significantly higher temperatures. 

Several studies have examined the case for a FlP-effect in 
other cool stars. Laming, Drake & Widing (1996) looked for a 
FlP-effect in the corona of e Eri. They found some evidence for 
a solar-like FlP-effect, but were unable to make a definitive detec- 
tion. More recently, Drake, Laming & Widing (1997) and Laming 
& Drake (1999) have identified solar-like FlP-effects in the coro- 
nae of a Cen and £ Boo. It has also been found, however, that there 
is no apparent FlP-effect in the corona of Procyon (Drake, Laming 
& Widing 1995). 

In summary, to date, there is no compelling evidence for the 
operation of a FlP-effect in the chromospheres/transition regions 
of the Sun or any other late-type stars. If there is such an effect in 



the coronae of stars, however, it is difficult to understand why it 
would not occur in the lower parts of the atmosphere. The C Il/Si II 
discrepancy in the modelling of e Eri could be interpreted as evi- 
dence for a FIP enhancement of the Si abundance in the low atmo- 
sphere of up to a factor of 2.4, and an atmospheric model based on 
that assumption is constructed in Section 6 (model C). There are, 
however, significant uncertainties and several shortcomings for this 
explanation which are also discussed in Section 6. 



6 MODELS B AND C 

In this section, the discussion of model A (Section 5.3) is used 
to develop two more models (B and C). These models were con- 
structed in the same manner as model A (see Section 5.1) but are 
based on different EMDs for the transition region and include a 
small amount of macroturbulence. The differences between mod- 
els B and C relate to the possible interpretations of the C Il/Si II 
discrepancy identified in model A (see below). The complete at- 
mospheric models (A, B and C) are tabulated in Appendix B. 



6.1 Modifications to the turbulent velocity field 

For model A it was assumed that all the observed non-thermal mo- 
tion was due to microturbulence. Since the optically thick chromo- 
spheric line profiles are sensitive to the distinction between macro- 
turbulence and microturbulence, the radiative transfer calculations 
for model A can be used to improve the prescription for the turbu- 
lence in the new models B and C. 

In Fig. 6, the computed Mg II h and k lines are plotted with 
the observed line profiles from STIS. Clearly the central reversal in 
the model A profiles is too wide, indicating that the assumed mi- 
croturbulence is too large. To compensate for this, macroturbulence 
is invoked to explain some of the observed line width: by trial and 
error, it was found that a value of £ mllC ro = 8.5 km s _1 leads to 
much better profiles for the Mg II lines. Unfortunately, there are 
no suitable lines to derive the macroturbulence in other parts of the 
atmosphere and so it is simply assumed that the macroturbulence 
does not vary through the upper chromosphere and transition re- 
gion: thus for models B and C, the microturbulent velocity field 
used above \ogT e [K] — 3.8 is 

Cmacro (4) 

where £(T e ) is the non-thermal velocity specified by equation (3). 
The adopted value for the macroturbulence has little significance 
above the low transition region, where it leads to less than 10 per 
cent difference in the adopted microturbulence. 



6.2 Modifications to the EMD 

For both models B and C the EMD has been slightly modified at 
high temperatures to improve the agreement between some of the 
predicted and observed line fluxes. In addition, more significant 
modifications have been made at low temperatures; these modifi- 
cations are discussed further below. The complete EMDs for the 
models are compared with that for model A in Fig. 7, and the 
atmospheric models obtained from the EMDs are compared with 
model A and the earlier models by Thatcher et al. (1991) and Kelch 
(1978) in Fig. 8. 
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Figure 6. Comparison of computed and observed Mg II profiles. In each figure the solid histogram shows the data, the solid curve shows the model A 
calculations, the dashed curve shows the model B calculations and the dotted curve shows the model C calculations. The narrow dips in the observed profiles 
are due to interstellar absorption. 
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Figure 7. Comparison of mean emission measure distributions. The solid curve is for model A, the dashed curve is for model B and the dotted curve is for 
model C. In each case the mean emission measure Em 3 (T e ) is plotted. 




Figure 8. Comparison of atmospheric models (column mass density versus electron temperature). The left panel shows the complete models and the right 
panel shows the transition region in detail. The heavy dashed line is the Kelch (1978) model, the heavy solid line is the Thatcher et al. (1991) model, the light 
solid line is model A, the light dashed line is model B and the dotted line is model C. 
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6.2.1 Model B 

Model B has been constructed based on the assumption that the dis- 
crepancy between the computed and observed C and Si fluxes is the 
result of errors in the computed C II resonance line surface fluxes, 
primarily due to the simplistic one-dimensional geometry used here 
(see Section 5.3.3). For this model, the DEMD was significantly 
increased around \ogT e [K] — 4.0 and slightly reduced at higher 
temperatures. These alterations increase the computed Si II fluxes 
and slightly reduce the C II and Si III fluxes relative to model A. 

6.2.2 Model C 

Model C has been developed to investigate the alternative explana- 
tion of the C Il/Si II discrepancy in model A, that of a (possibly FIP- 
related) change in the Si abundance. For this model, the assumed Si 
abundance was increased by a factor of 2.4. If this increase is due to 
a FlP-related effect then one would expect that the abundances of 
other low-FIP elements should be increased by a similar factor. Of 
the other element considered here this affects Mg and Al - hence 
the Mg and Al abundances have also been raised by a factor of 2.4 
relative to those in model A. 

The DEMD for model C was reduced throughout the region 
\ogT e [K] < 4.5 (compared with DEMD-A) to allow for the 
change in the emission measure loci resulting from the higher low- 
FlP-element abundances. 

6.3 H Lyman line profiles 

Since the cores of the strong H I Lyman lines are heavily affected 
by interstellar absorption, the diagnostic power of these lines for 
studying the stellar atmosphere lies in examining the line wings. 
Fig. 9 compares the Lyman a, f3 and 7 line profiles computed from 
models B and C with the observations. The model profiles have 
been corrected for interstellar absorption using the H I column den- 
sities and velocities derived for the e Eri line-of-sight by Dring et 
al. (1997). 

For Lyman (3 and 7, models B and C both predict the profiles 
very well. Model B tends to slightly overpredict the line strength 
but only by 10 - 20 per cent. For Lyman a the agreement is less 
good - both models predict profiles with significantly weaker line 
wings than observed. This suggests that the models underestimate 
the importance of photon scattering from the line core to the wings. 
Since this effect is not evident in the Lyman j3 or 7 profiles (nor 
the Mg II profiles) it seems unlikely that this is solely the result 
of errors in the atmospheric velocity field - most likely it indi- 
cates the coherence fraction used in the PRD treatment of this line 
is too high. A detailed investigation of the Lyman a profile goes 
beyond the scope of this work, but to illustrate the influence of a 
smaller coherence fraction, Fig. 9 also shows the Lyman a profile 
computed from model B adopting a frequency-independent, depth- 
independent coherence fraction of 7 = 0.9 (this is in contrast to 
the other calculations which use the physical 7 implemented in 
the Hubeny & Lites 1995 code). As expected, the constant-7 cal- 
culation predicts stronger line wings bringing it into much better 
agreement with the data. However, there are still some significant 
differences in the line shape which suggest that more redistribu- 
tion in needed in the line core but less to the far wings. The il- 
lustrative 7 adopted is rather low compared to theoretical calcula- 
tions and lower than the preferred Lyman a maximum coherence- 
fraction suggested by Vernazza et al. (1981) for their solar models 
(7 = 0.98). However, when comparing with the Vernazza et al. 



Table 10. Comparison of calculated fluxes (F c ) and observed fluxes (F ) 
for emission measure distribution models B and C. Following the model B 
fluxes, values are given for computations from a model identical to model 
B saving that the Asplund et al. (2004) relative abundances are adopted. F c 
given in units of 10~ 13 ergs cm~ 2 s — 1 . The top part of the table gives the 
density-insensitive lines, the middle gives the density-sensitive line and the 
lower part gives the FUSE lines. 



Ion 


A (A) 


F c 




(F c - F )/F x 100% 






Model B 


Model C 


Model B 


Model C 


N V 


1240 


1.27 (0.92) 


1.27 


-10 (-34) 


-10 


OV 


1371 


0.078 (0.053) 


0.078 


-25 (-50) 


-25 


Aim 


1854 


0.204(0.162) 


0.218 


-60 (-68) 


-57 


So 


1255 


0.190 (0.123) 


0.100 


-2 (-37) 


-48 


S iv 


1406 


0.090 (0.058) 


0.088 


21 (-21) 


19 


NIV 


1486 


0.069 (0.050) 


0.068 


52(11) 


52 


O in 


1663 


0.213 (0.144) 


0.196 


19 (-17) 


9 


OIV 


1400 


0.815 (0.551) 


0.814 


129 (55) 


129 


OV 


1218 


0.671 (0.454) 


0.671 


-1 (-33) 


-1 


NlII 


991 


0.452 (0.327) 


0.409 


37 (-1.0) 


24 


O VI a 


1035 


3.94 (2.66) 


3.94 


-43 (-61) 


-43 


S IV 


1070 


0.240(0.155) 


0.238 


156 (65) 


153 


S VI 


933 


0.100(0.065) 


0.100 


-58 (-73) 


-58 



a Flux likely to be underestimated since the line forms partially outside the 
temperature range covered by the models. 



(1981) value, it is noted that 7 is expected to be lower in e Eri than 
the Sun since the density is higher. 



6.4 Computed fluxes 

The new models (B and C) have been used to compute fluxes for all 
the lines considered for model A. Table 8 gives the fluxes computed 
from radiative transfer calculations and Table 10 gives the fluxes for 
the high temperature multiplets for which radiative transfer calcu- 
lations have not been performed (the model fluxes for these lines 
are obtained directly from the DEMD on the assumption that the 
multiplets are effectively optically thin). In both tables, the model 
B fluxes are followed (in parenthesis) by fluxes computed from a 
model identical to model B saving that it adopts the Asplund et al. 
(2004) relative abundances rather than those of Grevesse & Sauval 
(1998). 

The observed CI I 1351-A line provides a useful probe of the 
C II 1335- A radiation field within the atmosphere since the CI I line 
forms by fluorescence, pumped by the 1335- A transition (Shine 
1983). A very simple three level CI I model has been used to com- 
pute the 1351-A flux. The model contains only the ground con- 
figuration 2 P° states and the excited 3p 4 ( 3 P)4s 2 P 1/2 level. The 
only physical processes included in the model are the two radiative 
transitions P3/2 - P1/2 (the transition pumped by C II) and the 
observed transition ( 2 P°/ 2 ~~ 2 Pi/2)- The oscillator strengths were 
taken from Kurucz & Bell (1995). For the calculation the C II ra- 
diation field was taken from MULTI calculations with the full C 
model (Section 5.2.4). Table 11 gives the computed CI I flux for 
both models B and C. It is to be expected that the CI I flux may 
be overestimated since no levels of ionised CI I are included in the 
calculation (i.e. it is assumed that the CI I ion fraction is 1.0). 
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Figure 9. Comparison of computed and observed H I profiles. In each figure the solid histogram shows the data, the dashed curve shows the model B 
calculations, the dotted curve shows the model C calculations and the long-dashed curve shows the model D calculations. In the left panel (Ly a), the 
triple-dot-dash curve shows the model B result if a frequency-independent, depth-independent coherence fraction of 0.9 is adopted (see text). 



Table 11. Comparison of calculated fluxes (F c ) and observed fluxes (F ) 
for CI I in models B and C. F c given in units of 10~ 13 ergs cm~ 2 s _1 . 



Ion 


A (A) 


F c 


(F c - F )/F x 100% 






Model B Model C 


Model B Model C 


CI I 


1351.656 


0.231 0.110 


16.7 -80.0 


6.5 


Discussion 







The computed fluxes given in Tables 8, 10 and 1 1 are now used to 
assess the successes, failures and relative merits of models B and 
C. The following points are of particular note: 

(i) In general, both models B and C predict fluxes for the high 
temperature lines (Si III and hotter) that agree well with the obser- 
vations, given the combined uncertainties in the measured fluxes 
and atomic data (in particular the ionisation balance). 

(ii) The Mg II lines are predicted well by both models, both 
in terms of line profile and flux. The flux ratio is slightly differ- 
ent from that observed, possibly indicating that there is insufficient 
opacity in the models. It should be noted that if inhomogeneity af- 
fects the C II lines (as postulated for model B), the Mg II lines may 
also be affected - but without multi-dimensional calculations the 
scale and nature of such an effect is unknown. 

(iii) The Si II 1810-A multiplet is well predicted (within 20 per 
cent) by both models. This is the most important of the Si II mul- 
tiplets since it is the least sensitive to the absence of collision rates 
between the 4 P and more highly excited energy levels. 

(iv) The higher excitation Si II multiplets are predicted much 
less accurately. However, as discussed in Section 5.3.1 these may 
be significantly affected by inadequacies in the Si II atomic data 
which may account for differences of 50 per cent or more. 

(v) The C II resonance line fluxes (1335 A and 1037 A) are sig- 
nificantly overestimated by model B (~ a factor of two) but pre- 
dicted quite well by model C. As discussed in Section 5.3.3 the 
overprediction may be the result of the simplistic geometry as- 
sumed in the models. 

(vi) Although the uncertainties are large, the C II 1010-A mul- 
tiplet fluxes are not significantly overpredicted by model B. They 
are slightly underpredicted by model C. 

(vii) The CI I fluorescent flux is significantly underpredicted for 
model C. Model B predicts this line flux well. 

(viii) The S II flux is more accurately predicted by model B than 

C. 

(ix) If the Asplund et al. (2004) relative abundances are adopted, 
the model C II resonance line fluxes are closer to those observed 
but are still overpredicted. Adopting the new abundances does not 



significantly improve the agreement amongst the high-temperature 
lines since the relative abundances of C to N and O have not 
changed very much from those given by Grevesse & Sauval (1998). 

As discussed previously, most of the exceptions to (i) can 
be explained by inadequacies in either the ionisation balance for 
alkali-like ions (Al III, S VI) or the atomic data. However, the ten- 
dency of the density-sensitive lines to be overpredicted may also 
indicate the presence of density inhomogeneities and, therefore, the 
need for multi-component modelling (see Section 4.3). 

Points (vi), (vii) and (viii) all argue in favour of model B over 
model C. The only significant advantage of model C over model B 
is that it predicts the C II resonance line surface fluxes more accu- 
rately. However, model C has one more free parameter than model 
B (the FlP-enhancement factor) and so it is not surprising that bet- 
ter agreement can be found. Point (vii) is particularly crippling for 
model C since it suggests that the model does not reproduce the 
C II radiation field within the atmosphere. 

Thus it is concluded that model B is the better single- 
component model for the upper transition region of e Eri, and 
that its main limitation is the assumption of homogeneity. In fu- 
ture work inhomogeneous models and multi-dimensional radiative 
transfer must be employed to resolve the issues relating to scatter- 
ing in the C II resonance lines. 

The interpretation of the apparent discrepancy in the C Il/Si 
flux ratios in terms of a variation in the abundance ratio (model C) 
is judged to be inferior. Model C is, however, a useful constraint 
on FlP-related abundance enhancement in the transition region of 
e Eri: any enhancement in the Si abundance is constrained to be no 
larger than that adopted in model C (a factor of 2.4). 



7 A SIMPLE MULTI-COMPONENT ATMOSPHERIC 
MODEL 

It is well known from direct observations of the Sun and studies of 
other stars (e.g. the K giant a Tau, McMurry & Jordan 2000) that 
the atmospheres of cool stars are inhomogeneous. As mentioned 
in Section 5.3.3, Sim & Jordan (2003b) concluded that emitting 
material occupied only a limited fraction (~ 20 per cent in the 
mid-transition region) of the available volume in the atmosphere 
of e Eri. All the models discussed above (A, B and C) are single- 
component and assume a filling factor of 100 per cent. One im- 
portant consequence of a small filling factor is that line-of-sight 
opacities will be higher since the plasma is confined by a smaller 
cross-sectional area. To fully investigate multi-component atmo- 
spheric models requires multi-dimensional radiative transfer which 
goes beyond the scope of this paper, however, in this section a crude 
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two-component model will be presented and contrasted with the 
homogeneous models discussed above. 

For the two-component model (model D) it is assumed that 
the atmosphere has a uniform filling factor of emitting material, 
fixed at 20 per cent (i.e. A(r)/A,(r) = 0.2). The apparent EMD 
of the emitting material is the same as that used for model B. For 
model D the intrinsic EMD is therefore larger by a factor of 5 and 
hence leads to larger line opacities. The remainder of the volume is 
assumed to be filled by a completely dark second component. One- 
dimensional radiative transfer calculations are performed for the 
emitting component of the model neglecting the dark component. 
This model can only be regarded as exploratory since it neglects the 
probable process of photon leakage from the emitting component 
to the dark component but it is useful to investigate the effects of 
higher line opacities. 

Fluxes computed from radiative transfer calculations with 
model D are given in Table 8. Since the EMD for model D is iden- 
tical to that of model B, the model D fluxes for lines assumed to 
be optically thin (Table 10) are the same as those for model B. The 
CI I flux for model D is also virtually identical to that of model B 
(Table 11). 

Comparing the model B and D fluxes in Table 8 reveals sev- 
eral interesting changes, mostly in resonance lines that form at rel- 
atively low temperatures (these being the lines for which opacity 
effects are most important). The most dramatic differences are in 
the absolute fluxes of the strongest lines, those of H I and Mg II. 
Increased opacity reduces the fraction of emitted photons which 
eventually escape the atmosphere in these lines. The effect is great- 
est in the H Lyman lines where the model D fluxes are reduced 
by a factor ~ 2, relative to those from model B (see Fig. 9). Al- 
though this shows that the 1-D models presented here (B and C) 
predict the Lyman lines more accurately than the two-component 
model (D), it must be noted that the very simple implementation 
of model D neglects photon escape from the emitting component 
to the second component. Therefore it is expected that model D 
will underpredict optically thick line fluxes relative to a more so- 
phisticated, multi-dimensional radiative transfer calculation for a 
multi-component atmosphere. 

The Mg II lines are less affected than those of hydrogen but 
still decrease significantly in strength (by a factor of ~ 1.6). The 
computed ratio of the two Mg II lines changes in the sense which 
brings it closer to the observed ratio, but the change is very small. 

For lines of all other ions, the changes in absolute flux are 
quite small, at most ~ 10-20 per cent. For some pairs of Si II 
lines with common upper levels (i.e. the 1526.7/1533.4-A pair and, 
in particular, the 1808. 0/1817. 5-A pair) the line ratios from model 
D are in better agreement with the data that those from model B. 
The provides some evidence that the line-of-sight opacities sug- 
gested by the two-component model are more realistic than those 
in the homogeneous models. However, the 1309.3/1304.4-A flux 
ratio from model D agrees less well than that from model B (note, 
however, that the absolute fluxes for these lines are generally in 
poorer agreement with the data than those of the other lines consid- 
ered here). 

The C II 1334.7/1334.5-A flux ratio changes from 2.55 in 
model B to 2.71 in model D, both these values being consistent 
with the observations, while the 1036. 3/1037. 0-A ratio is in sub- 
stantial disagreement with the data for both models. 

It is concluded that the simple two-component model pre- 
sented here is not significantly superior to the homogeneous mod- 
els discussed above, although there are some improvements for the 
modelling of the most important Si II line ratios. A more complete 



treatment incorporating multi-dimensional radiative transfer effects 
is needed to fully investigate a multi-component model; this is de- 
ferred to later work. 



8 SUMMARY AND CONCLUSIONS 

Using measured ultraviolet emission line fluxes and widths, new 
single-component models for the upper chromosphere/transition 
region of e Eri have been constructed; a preliminary model (model 
A) and two further models (models B and C) which address some 
of the inadequacies of the first model. To date, these models 
are the most detailed that have been constructed for the chromo- 
sphere/transition region of a late-type main-sequence star, other 
than the Sun. A simple two-component model (model D) has also 
been presented. The predictions of that model were only subtly dif- 
ferent from those of the equivalent single-component model, sug- 
gesting that the 1-D modelling does provide useful average models. 

For the emission measure analysis involved in the construc- 
tion of the models new ionisation balance calculations have been 
performed which attempt to account for the density sensitivity of 
dielectronic recombination rates. Although not definitive, these cal- 
culations highlight the extreme sensitivity of the analysis of certain 
ions (generally the alkali-like ions) to errors in the recombination 
rates and illustrate the need for the publication of reliable recombi- 
nation rates at finite density. 

In most respects, line fluxes/profiles computed with the 
favoured model (model B) are in agreement with the observations, 
allowing for modest measurement uncertainties and reasonable er- 
rors in the adopted atomic data and ionisation/recombination rates. 
The most significant failures of the model are its inability to com- 
pute reliable fluxes for the optically thick resonance lines of C II 
and the spin-forbidden transition region lines of several ions (most 
notably O IV). Barring large errors in the atomic data, the most 
probable explanation for the discrepant transition region forbidden 
line fluxes lies with density variations in the atmosphere. The C II 
resonance lines have been examined in detail and used to place 
a constraint on possible FlP-related abundance variations within 
the upper atmosphere. However, it has been concluded that the 
anomaly associated with these lines is most likely to be related to 
scattering in an inhomogeneous atmosphere. 

Thus the models presented here represent the limit to what 
can be achieved with single-component models of the outer at- 
mospheres of late-type stars, and the single-dimensional radiative 
transfer calculations associated therewith; such models can account 
for most of the observed properties of ultraviolet emission lines, 
but modern spectra are of sufficient quality to indicate some inad- 
equacies. Evidence of inhomogeneity comes as no surprise; direct 
observation of the Sun shows us that stellar atmospheres contain 
structure. However, at present we have no direct means of observ- 
ing inhomogeneity in other main-sequence stars and so it can only 
be understood by modelling spectra. 

Theoretical studies involving two-component models are al- 
ready underway (e.g. Cuntz et al. 1999). In the future con- 
straints on modelling of this type must be provided by semi- 
empirical multi-component modelling of spectra involving multi- 
dimensional radiative transfer calculations (codes capable of such 
multi-dimensional calculations are now available e.g. Uitenbroek 
2001). 

The indications of inhomogeneity in the upper chromo- 
sphere/lower transition region of e Eri dovetail with our analysis 
of the upper transition region/corona of the same star (Sim & Jor- 
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dan 2003b), in which the filling factor of emitting material in the 
mid-transition region was found to be ~ 20 per cent. Given the 
current absence of a definitive understanding of the heating mech- 
anisms which operate in the lower transition region and below, a 
similar argument cannot be applied to directly estimate the filling 
factor in these parts of the atmosphere, but if the filling factor is as- 
sociated with magnetic fields (as expected) then solar observations 
would suggest that the filling factor should become smaller deeper 
in the atmosphere. However, a full investigation of this must await 
quantitative 3-D modelling of the optically thick lines which can 
provide information about the inhomogeneity; it is to be hoped that 
in the near future such modelling can be combined with advances 
in theoretical studies to lead us to a clearer understanding of the 
structure and heating of the outer atmospheres of late-type stars. 
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Al. Logarithmic ionization fractions computed for C (see Section 3). 



log T e (K) Ionization fraction 
C II C III C IV 
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-5.29 


-2.72 
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-3.60 
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-7.33 
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-4.89 


-2.56 



APPENDIX A 

Tables Al - A4 give the computed logarithmic ionization fractions 
for C, N O and Si, respectively. See Section 3 for details of these 
calculation. 



APPENDIX B 

Tables Bl - B3 give parameters for the atmospheric models A, B 
and C (as presented in Sections 5 and 6, respectively). 
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A2 Logarithmic ionization fractions computed for N (see Section 3). 



logT e (K) Ionization fraction 
N III N IV Nv 
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A3 Logarithmic ionization fractions computed for (see Section 3). 



log T e (K) Ionization fraction 

OlII OIV Ov OVI 
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A4 Logarithmic ionization fractions computed for Si (see Section 3). 
log T e (K) Ionization fraction 
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Bl Model A. m co ; is the column mass density, T e is the electron tem- 
perature, N e is the electron number density, Nh is the hydrogen number 
density and £ m j cro is the most-probable microturbulent velocity. 



log m co i 
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B2 Model B. m coi is the column mass density, T e is the electron tem- 
perature, N e is the electron number density, Nu is the hydrogen number 
density and £ m i cro is the most-probable microturbulent velocity. 



log m col 


logT e 


logiV e 


logN H 
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1 1.41 
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10.28 


-4.50164 
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11.56 
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-4.50037 


4.110 


11.54 


11.61 


9.14 


-4.49767 


4.070 


11.57 


11.67 


8.55 


-4.49026 


4.030 


11.59 


1 1.73 


7.93 


-4.47013 


3.990 


11.60 


11.83 


7.29 


-4.43521 


3.950 


11.55 


11.95 


6.61 


-4.37097 


3.910 


11.45 


12.11 


5.89 
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3.870 


1 1.29 
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5.10 
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4.40 
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1.34 
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15.51 


1.27 
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15.68 


1.25 


-1.06000 
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11.60 


15.89 


1.25 
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3.619 
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16.08 


1.25 
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16.33 


1.25 
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3.619 


12.08 


16.55 


1.25 
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3.621 


12.26 


16.79 


1.25 


0.06000 


3.626 


12.43 


16.99 


1.25 


0.26000 


3.641 


12.66 


17.18 


1.25 


0.43000 


3.666 


12.93 


17.33 


1.25 


0.59893 


3.704 


13.24 
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3.736 
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17.54 
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0.79292 
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0.87069 


3.788 


13.92 


17.65 


1.25 



B3 Model C. m co ; is the column mass density, T e is the electron tem- 
perature, N e is the electron number density, Nh is the hydrogen number 
density and £ m j cro is the most-probable microturbulent velocity. 



log m co i 
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logTe 
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-4.54760 
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10.27 
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-4.54695 


5.090 


10.64 


10.58 


19.15 


-4.54685 


5.01 1 
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A C A f A O 

-4.54648 


4.649 


11.05 


10.98 


16.16 


A C A if A A 

-4.54644 


4.61 1 


11.08 


1 1.02 
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14.20 
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8.24 
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